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1. Introduction 


Recent data from Au + Au reactions at yj snn — 130 and 200 GeV from RHIC show 
striking difference between baryons and mesons in the intermediate transverse momentum 
region 2 < p± < 5 GeV. First, elliptic flow scales with constituent quark number [[T1I21 El IH 
The flow per constituent quark as a function of transverse momentum per constituent 
quark, v% adron ( pt/u)/ n (n — 2 for mesons, 3 for baryons), is a universal function for all 
hadron species (within experimental uncertainties). Second, nuclear suppression is weaker 
(Raa is larger) for baryons than for mesons[ HI [21 0 EL Si El 21- The most promising 
mechanism proposed to explain both phenomena is parton coalescence [ HIH HT11121 131 1T4I 

fToi m on . 

In the coalescence model, mesons form from a quark and antiquark, while baryons from 
three quarks or three antiquarks (contributions from higher Fock states are small [ 1TH1). 
Most versions of the model are based on variants of the simple “coalescence formula” [ 
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that gives the hadron spectra in terms of the constituent phasespace distributions /) on 
a 3D spacetime hypersurface and the Wigner-transform of hadron wave functions W. 
(Arbitrary 3D hypersurfaces[ [HE 20J and quark correlations[ EH are straightforward to 
accommodate.) This simple approach can reproduce the particle spectra at RHIC quite 
well, with coalescence hadronization from a thermalized quark-antiquark plasma] 13] H3jj 
and an additive fragmentation contribution of quenched high-p^ jets[ 22]- The scaling of 
elliptic flow with constituent quark number was also explained in this framework [113. 131 . 





2 


Denes Molnar 


Nevertheless, several important questions are still open. At low pr, the coalescence for¬ 
mula violates [|2l El unitarity. The yield in a given coalescence channel scales quadrat- 
ically/cubically with constituent number, moreover, the same constituent contributes to 
several channels (including fragmentation in certain schemes). Energy conservation is 
only approximate because of the on-shell treatment and neglect of binding energies. Also, 
though the final entropy is comparable to the initial one (due to hadron decays) [ |E], it 
is problematic that entropy does decrease temporarily during the coalescence process. 

Most importantly, the simple parametrizations assumed for the quark phase space distri¬ 
butions are inconsistent with dynamical models, such as parton cascades or hydrodynam¬ 
ics. Spatial inhomogeneities and dynamical phasespace correlations distort coalescence 
predictions in a crucial way. In particular, the observed quark number scaling becomes 
highly nontrivial [IM|, and a large baryon/meson ratio enhancement is difficult to achieve[ 

Ej. 

Below I review and illustrate the difficulties, using one of the main dynamical ap¬ 
proaches, covariant parton transport theory. 


2. Covariant parton transport theory 


Covariant parton transport theory is the incoherent, particle (short-wavelength) limit of 
quantum-chromodynamics. Its main advantage is that it is applicable out of equilibrium 
and models freezeout self-consistently. However, it cannot describe phase transitions 
(without coupling to classical fields). 

One of the simplest but nonlinear forms of Lorentz-covariant Boltzmann transport 
theory involves on-shell phase space densities {fi(x,p}} that evolve with elastic 2 —» 2 [ 
EE El EH El El and inelastic 2 —* 2 [ El rates as 
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| Ad | 2 is the polarization averaged scattering matrix element squared, the integrals are 
shorthands for f = f d 3 p a /(2E a ), gi is the number of internal degrees of freedom for 

a 

species i, while f a ^ = fi(x,p a ). The source terms {S)(x,p)} specify the initial conditions. 

Eq. (0) could in principle be extended for bosons and/or for inelastic processes, such 
as gg ggg [EElilSJ- However, with the new nonlinearities these extensions introduce, 
it is very challenging to maintain Lorentz covariance numerically at opacities expected at 
RHIC. Off-shell variants of parton transport also exist, such as the VNI/b model[ [Ej. 

Based on the local mean free path A (x,s) = l/[n(:r)(7(s)], transport theory naturally 
interpolates between ideal hydrodynamics (A —> 0) and free streaming (A —> oo). The 
most relevant quantity that characterizes the whole evolution is the transport opacity 
X = f dzp(z)a tr [ EJ) which is the average number of collisions per parton multiplied by 
the ratio of the transport and total cross sections (the latter is the efficiency of momentum 
transfer in a single collision). In a near-equilibrium situation, these parameters can be 
related to transport coefficients, such as the shear viscosity or diffusion constants. 

Transport theory provides important information about the properties of the partonic 
medium created at RHIC. The large elliptic flow observed[ HUES EH] indicates a strongly- 
interacting, opaque parton system[ IE] , with dN g /dg(b = 0) x <j gg ~ 1000 x 45 mb, about 
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Figure 1. Results for Au + Au at y/s = 200 A GeV at RHIC with 6 = 8 frn from MPC[ 
HO!. a) Transport opacity dependence of elliptic flow as a function of px from covariant 
transport theory [EE!; b) comparison of elliptic flow as a function of px from covariant 
transport theory and ideal hydrodynamics [EH]- 


15 times above the perturbative estimate (see Fig. |T}a) . These conditions are still different 
from an ideal fluid because dissipation reduces elliptic flow by 30 — 50% relative to the 
ideal fluid limit [EH] (Fig. Eb). At such high opacities, the shear viscosity of the parton 
plasma p ~ sAT/5 [ S3 (estimated from kinetic theory) is indeed very small, close to 
the conjectured lower bound p min = s/(47t) [ SEJ (s is the entropy density). However, in 
heavy ion collisions, the small viscosity is compensated by large gradients, resulting in 
significant dissipative effects. Estimates based on viscous (Navier-Stokes) hydrodynamics 
also support this conclusion[ 5H1 . 

Though even the large opacities at RHIC are insufficient to equilibrate heavy quarks, a 
large charm quark elliptic flow is still expected above px ~ 2 — 3 GeV (Fig. El) from several 
studies (covariant transport predictions from MPC[J3II! 4DJ, reinforced by calculations from 
AMPT[ 3Tj and also a Fokker-Planck approach[132]). It is very exciting that preliminary 
indirect data from RHIC do indicate a sizeable charm elliptic flow[ ITT! 133] . These data 
are on the elliptic flow of electrons from D meson decays 1%*) —► Kue, which based on 
model studies follows closely the D meson elliptic flow[135, 13T] (see Fig. Et>)- 

At high px, covariant transport shares a lot of similarity with parton energy loss models [ 
E2]> except that 2^2 transport gives incoherent, elastic energy loss. However, at the 
high opacities at RHIC the dynamics is much richer. Not only energy loss, but significant 
energy gain is also possible in multiple collisions [ |3B]. This “plasma push” process, a 
shadow of uear-hydrodynamic behavior at low px , could play a role even at px ~ 10 GeV, 
and is the reason why elliptic flow from the transport decreases very slowly at high px- 
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Pt [GeV] 




Figure 2 . a) Top panel: parton species dependence of elliptic flow in Au + Au at y / Jnn = 
200 GeV at RHIC as a function of pt, for b = 8 frn and ~ 7 times perturbative opacities, 
from the covariant MPC model [EDI; bottom panel: preliminary non-photonic electron 
V2 (pt) data from RHIC [ 143j . b) Top panels: parton u 2 (jpr) for minimum-bias Au + Au 
collisions at yj s^n — 200 GeV from the AMPT model; bottom panels: comparison of D 
meson w 2 and electron u 2 from D decays [ EfP . 


3. Dynamical effects and elliptic flow scaling from coalescence 

To illustrate the effect of the dynamics on elliptic flow scaling from coalescence, it 
is useful to first ignore variations of phasespace distributions on length and momentum 
scales corresponding to a typical hadron (~ 1 fin and ~ 200 MeV). In this case[ ’ 21] . 
W ~ 5 3 (Axij)S 3 (Apij) and the phasespace distributions of mesons and baryons from 
coalescence ( a/3 —> M, a/3 7 —»• B) are 

/ b ( x , p ) = —- 9B f a (x,p/3) fp(x,p/3)f 1 (x,p/3) 

9a 9/3 9~1 

Im(x,p) = ( 2?r ) 9m f a (x,p/2) fp(x,p/2) . (3) 

9a 9(3 

g is the degeneracy of the particle (spin and color), and I considered constituent quarks 
of (roughly) equal mass, which share the total hadron momentum (roughly) equally. (I 11 
hadrons that also contain a heavy quark, e.g., D mesons or the A c , the heavy quark carries 
most of the momentum [tT5].) 
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Assume for simplicity that all quarks have identical phase space distributions 


f q {x,p T ,y = 0) = n q (x,p T ) 


1 + 51 2v n,q(x,p T ) cos(n 0 ) 


n =1 

where n q is the local density of quarks with transverse momentum magnitude px- 
small local anisotropies 


(4) 

For 


V2,q\ < 1 , \Vk,qVl,q\ < K q\ , 


(5) 


the usual assumption to obtain elliptic flow scaling, the hadron elliptic flows from © are 


v 2 ,m(Pt) ~ 
V2 ,b{Pt) ~ 


2 ( n 2 q (x,p T /2 ) v 2:q (x,p T /2)) x 

(n 2 q (x,p T /2)) x 

3 (n q (x,p T / 3) v 2 , q {x,p T /3)) x 

( n q( X ,p T /3))x 


( 6 ) 


with (A(x)) x = f d 3 xA(x) (spatial average). Though each small spatial region with its 
own local quark elliptic flow contributes a local hadron elliptic flow that scales with quark 
number, v 2 ,had( x iPT) ~ nv 2>q (x,px/n), quark number scaling does not follow in general 
for the measured (spatially averaged) elliptic flow, i.e., v 2t had{PT ) ^ nv 2>q (px/n). This is 
because constituent phase space densities appear nonlinearly, raised to second and third 
power. 

Constituent quark scaling arises only for special classes of quark phase space distribu¬ 
tions. Two classes where the spatial dependence in © cancels out trivially are i) when 
the local density is spatially uniform n(x,px) = n(px) (as assumed in [jT^J); and ii) when 
the local flow anisotropy is spatially uniform v 2 (x,px) = v 2 (px) (the assumption in [ 113] ). 

Unfortunately, neither of the two cases is realistic. Dynamical models yield spatially 
nonuniform local densities and nonuniform flow anisotropies. Local density variations have 
been studied in connection with the pion interferometry puzzle from hydro dynamics [ I47| 
and several transport models[ 14811491 5U] (Fig. Eh is an example from covariant parton 
transport). Local flow anisotropies are also spatially nonuniform as can be seen in Fig. Eb 
that shows (momentum) azimuthal angle <fi distributions at midrapidity for 1 < px < 2 
GeV, for Au + Au at top RHIC energy, 6 = 8 fm, averaged over four wedges in the 
transverse coordinate plane ip x E [kir/8,(k + 1 ) 7 t/ 8 ], with k = 0,1,2,3 [El). Clearly, 
partons emitted from in-plane regions move predominantly in-plane and, therefore, have 
positive elliptic flow v 2 > 0; while partons coming from out-of-plane regions move in out- 
of-plane direction and have negative v 2 < 0. This is very similar to the hydrodynamic 
expectation[ 5JJ based on radially boosted fireballs. 

Fig. Eb also demonstrates that the condition of small local anisotropies (j5|) . which 
is necessary for general quark number scaling, is not satisfied either. Instead of small 
harmonic modulations over a uniform background, the distributions are strongly peaked 
because high-p^ particles can only escape from a surface layer of the reaction region. 
In this case, local anisotropies from coalescence follow a unique power-law scaling [El 
\ v had{pT) — \v 2t q(px/n) | x , instead of the linear one derived in m. 

Therefore, quark number scaling of elliptic flow requires a highly nontrivial, fortuitous 
interplay between variations in local parton density and nonlinear couplings between large 
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Figure 3. Results from MPC[0UJ for Au + Au at y ^snn = 200 GeV with b = 8 fm. a) 
Transverse positions of partons at freezeout as a function of pt, for final rapidities \y\ < 2 
(scatter plot). The ellipses guide the eye. b) 0 distributions at midrapidity for 1 < p T < 2 
GeV, averaged over four wedges in the transverse plane <p x E [k7i/8, (k + 1)7t/8], k = 0...3 

w- 


local flow coefficients. Several classes of phase space distributions have been explored 
recently [ 52 f53] and solutions have been found that scale approximately. 

One particularly simple example[ 53], with large local flow anisotropies, is when N par- 
tons move along the 0 = 0 axis and (1 — a)N partons move along 0 = 7r/2(a-Cl). If the 
local volumes for both components are identical, the hadron yields are Nhad, oc V ( N/V) n -\- 
V((l—a)N/V) n , and the elliptic flows do scale v 2 = [(N/V) n — ((1 — a)N/V) n ]/[(N/V) n + 
((1 — a)N/V) n } ~ na/2. However, it is easy to show that scaling crucially depends on 
the assumption of identical volumes. Different local volumes, V and V' , give hadron 
yields N had oc V(N/V) n + V'((l — a)N/V') n and elliptic flows that do not scale , v 2 = 
[( V/V'y- 1 -(1 -a) n ]/[( V/V'y- 1 + (I-a) n ] « [( V/V') n - 1 -1 + na)}/[(V/V') n ~ 1 + 1 +na)}. 

This underscores the importance of going beyond simple parameterizations and using re¬ 
alistic dynamical models to study parton coalescence in heavy-ion collisions. 

It is important to realize that quark number scaling from coalescence necessitates 
nonequilibrium dynamics[ 2Tj. For thermal constituent distributions, coalescence reduces 
to statistical hadronization and, therefore, momentum anisotropies in that case can only 
depend on particle mass, and not quark number. This is so even for nonequilibrium quark 
distributions that are of the form f(x,p) = g(p fl u IJ '(x),x). Therefore, besides the break¬ 
down of hydrodynamic behavior, quark number scaling indicates a marked departure from 
pure “flow-like” space-momentum correlations as well. 

4. Dynamical coalescence approach 

One way to study realistic dynamical effects is to combine nonequilibrium parton trans¬ 
port theory and coalescence [E7J. Such an approach requires overcoming the limitations of 
fT]) and proper book-keeping to ensure that each parton participates in one hadronization 
channel (no multiple counting). 

The simple coalescence formula m relies on the assumption that interactions between 
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quarks that are not in the same hadron cease suddenly on some 3D hypersurface. However, 
self-consistent freezeout from transport approaches gives diffuse 4D freezeout distribu¬ 
tions [HBJ US ESI that cannot be well approximated with a hypersurface. The coalescence 
formalism was extended to such a case by Gyulassy, Frankel and Render (GFR) in [ Ell 
(for weakly-bound states). Their result is the same as (fTj) . except that the weight W is eval¬ 
uated using the freezeout coordinates ( t , x ) of constituents. When taking Ax, the earlier 
particle needs to be propagated to the time of the later one, e.g., AT 12 = X\ — T 2 + (^ 2 — G)d 
if ti < t 2 . The origin of this correction is that a weak bound state survives only if none 
of its constituents have any further interactions. For baryons, the generalization involves 
propagation to the latest of the three freezeout times. 

To facilitate proper book-keeping, it is common in transport models[ ES to utilize 
box Wigner functions W = TLy Q(x m — Aiy,1)0 (p rn — Ap t j |). This way (HJ) has a 
simple interpretation: if IT = 1 (and the quantum numbers match) the hadron is formed, 
otherwise it is not (W = 0). If several coalescence final states exist for a given constituent, 
one is chosen randomly. Partons that do not find a coalescence partner would propagate 
“freely” and are therefore fragmented independently. 

Fig. 0 shows pion and proton elliptic flow at RHIC for Au + Au at sJsnn = 200 GeV 
with 6 = 8 fm and ~ 7 times perturbative opacities (a 99 ^ 99 = 10 mb), from the dynamical 
coalescence approach[ TZj based on the covariant transport model MPG[ET)|. The elliptic 
flow of direct pions and protons from coalescence (left panel) is somewhat smaller than 
the parton flow scaled by constituent quark number. The difference is due to genuine 
dynamical effects that are much larger than few-percent corrections nonlinear in u 2 or 
those due to higher-order flow anisotropies[II5j. Nevertheless, in view of the discussion in 
the previous Section, it is quite remarkable that the end result is only a modest ~ 20% 
and ~ 30% flow reduction for mesons and baryons. 

On the other hand, fragmentation contributions sharply reduce the hadron u 2 as shown 
in the right panel in Fig. El Unlike the flow amplification from coalescence [ JI3J, fragmen¬ 
tation smears out the anisotropy because hadrons from the parton shower are not fully 
collinear (nonzero jet width). Because only a smaller fraction of parton finds a coalescence 
partner, a large fraction fragments, including essentially all partons above pr > 2.5 GeV, 
as can be seen in Fig. Eh. Coalescence does become more and more important at low pt 
but it involves only about one in every three partons. 

Despite the large fragmentation component, pion and proton u 2 still scale approximately 
relative to each other - up to small differences of ~ 10%, perhaps as much as 30%, which 
could be tested against high-precision Run-4 data from RHIC (better theory accuracy in 
th epx/n > 1 — 1.5 GeV region is certainly desirable). However, quantitative interpretation 
of the data becomes more complicated because the scaled hadron flows underpredict the 
real parton u 2 (pt) by factor of 2 or more. 

Dynamical effects on the baryon/meson ratio are much more striking as can be seen in 
Fig. Eb Here the yield enhancement from coalescence is characterized by R coa h the ratio 
of the final spectra with hadronization via combined coalescence and fragmentation to 
that with hadronization via fragmentation only. Parton coalescence enhances both pion 
and proton yields, and hence Raa, by as much as a factor of three in the “coalescence 
window” [CEDES! i-5 < Pt < 4.5 GeV. The additional hadron yield comes dominantly 
from partons with 0.5 < Pt < 2 GeV (see Fig. Eb)- The enhancement, however, is 
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Figure 4. Quark number scaled elliptic flow from MPC[ j^Q] as a function of px for 
pions (open squares) and protons (filled circles) in Au + Au at yfs = 200 A GeV at 
RHIC with b = 8 frn and cr gg = 10 mb, with hadronization via combined coalescence and 
fragmentation (right), and for primary hadrons (without decays) from coalescence (left). 
The constituent n 2 (Pt) is also shown (solid lines). Taken from [ T71 . 


about the same for both protons and pions, and hence p/tt stays close to the value in 
p + p collisions. Even for a (rather unrealistic) scenario with immediate freezeout on the 
formation r = O.lfm/c hypersurface, which results in very high parton densities, the p/n 
ratio is enhanced relative to p + p only 1.5 — 1.7 times, much less than the factor fa 3 seen 
at RHIC at this centrality[Ej. 

5. Conclusions 

Parton coalescence is a promising approach to explain the striking baryon-meson differ¬ 
ence in the Au + Au data at intermediate 2 < pj < 5 GeV from RHIC. Studies based on 
simple parameterizations are quite successful but are inconsistent with dynamical models. 
A dynamical coalescence approach based on parton transport theory does confirm that 
coalescence is an important hadronization channel at RHIC. However, the dynamics sig¬ 
nificantly affects elliptic flow scaling and baryon/meson ratios. Until a reliable approach 
can successfully reproduce the data, “the jury is still out” regarding coalescence. 

Clearly, further studies are needed. The 2 —> 2 transport considered could be an 
oversimplification, or the formalism based on weakly-bound states may not apply well to 
the QCD coalescence process. Long-range correlations could make coalescence nonlocal 
in phase-space (unlike (fljl ). or help “optimally match” coalescence partners to maximize 
the yields. Precise experimental data on deviations from quark number scaling, flavor 
dependence of observables, and the transition to the pure fragmentation regime at high px 
will provide additional tests of coalescence models and further insights into the dynamics 
at intermediate px- 
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Figure 5. Results [ I7j from MPC for Au + Au at y/s = 200A GeV at RHIC with 6 = 8 
frn for o gg = 3 mb (dotted) and 10 mb (solid), or immediate freezeout at r = 0.1 frn/c 
(dashed-dotted line), a) Fraction of partons that fragment independently as a function of 
Pt ; b) pion and proton enhancement from coalescence as a function of p±. 
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